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Abstract. The phase-portrait of the second order differential equation: 



n 



is studied. Some results concerning existence, non-existence and uniqueness 
of limit cycles are presented. Among these, a generalization of the classical 
Massera uniqueness result is proved. 



1. Introduction 



In this paper we investigate the qualitative behavior of the differential equation: 



where P{x, y) — X]r=o fii^)y'' > ^^r some fixed n G N, which in the phase-plane can 
be rewritten as the differential autonomous system: 



Such systems have been proposed by H.I. Freedman and Y. Kuang ^ in the study 
of a Gause-type predator-prey model. 

If we assume some standard regularity assumptions on the functions fi{x), 
I G {0, . . . ,n}, for the uniqueness of Cauchy's problem, for instance let foix) — 



Al) g e Lip{R); 

A2) fieC{R), le{l,...,n}; 

and the hypothesis: 

B) xg{x) > for aU a; ^ 0; 

then the origin is the only singular point, and trajectories turn clockwise around it. 

The problem of existence and uniqueness of limit cycles for system will 
be considered. All examples and applications presented will involve polynomials, 
because this is, from our point of view, the more interesting situation; we never- 
theless observe that all our results hold in a more general setting than polynomial 
functions. 



(1.1) 



x + P{x,x) = 0, 



(1.2) 
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Clearly equation Hl.l|l is a particular case of the classical generalized Lienard 
equation: 

X + f{x, x)x + g{x) = , 
as well as the equivalent system: 

X 

y = -f{x,y)y - g{x) . 

is a generalization of system H1.2|l . 

We observe that a large number of results, in the huge literature concerning these 
generalized Lienard equations, are obtained with some assumptions on the growth 
of f{x,y), for instance: 

f{x,y) > ~M, 

being M some suitable positive constant. This because otherwise trajectories might 
be not continuable and the uniform boundedness of the trajectories is a crucial step 
in order to apply the classical Poincare-Bendixon Theory. 

This is not the case of system (|1.2|) where in general we will consider functions 
fi{x), Z £ {0, . . . , n} to be polynomials, hence most of the classical result cannot be 
used here. 

As far as we know system (|1.2(l has been studied only in two particular cases 
(the classical Lienard equation): 

Pi{x,y) = g{x) + fi{x)y, 

and 

P2{x,y) ^ g{x) + fi{x)y + f2{x)y'^ , 
largely studied by several authors in the last decade (see for instance |2 121 [SJ ^ ) . 

We believe that the main reason for this, is the existence of a well known trans- 
formation, discussed in sectional bringing equation: 

(1.3) X + h{x)x + f2{x)x^ + g{x) = , 

into some Lienard system. In this way all classical results of existence of limit 
cycles for the Lienard equation may be used. A second reason, for considering 
systems (|1.3(l . lies in the fact that its trajectories are always continuable. Of course 
this can be obtained via the transformation into a Lienard equation, but even in a 
direct way, considering the slope of the trajectories in the phase-plane. 

When one is interested in general cases, these properties are no longer valid. 

In this paper we will exploit some geometrical properties of the phase-portrait 
of system (|1.2|) and when P{x, y) = P2{x, y), we will not make use of the transfor- 
mation of the system into a Lienard one, but we will use a slight different system. 

The plan of the paper is the following. In section El we will study system H1.2II . 
splitting P into two parts: one exhibiting some suitable symmetries and consid- 
ering the other as a " perturbation" , then some properties will be enlighten using 
a comparison method. The main result of this section is the following Theorem, 
which gives a necessary condition for the existence of limit cycles. 

Theorem 1.1 (Non-existence). Let us consider the second order differential equa- 
tion: 

N M 

(1.4) i + Y, /2/„i(a;)i2'-i + ^ f2i{x)x^' + g{x) = , 

1=1 1=1 
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where N and M are positive integer. Assume regularity hypotheses A and B to he 
satisfied by the fi{x) 's and g{x), then systems ()1.4|l has not periodic orbits provided 
all the f2i-i{x) 's, I G {1, . . . , N}, never change sign. 

Section El deal with system (11.3(1 . As already mentioned the study is performed 
using a new differential system rather than transforming system into some 

Lienard equation. The main result of this section is the following: 

Theorem 1.2 (Existence). Assume hypotheses A and B to hold, let Fi{x) — 
lo fi{s)ds , G{x) — /q g{s) ds and assume: 

C) there exists S > s.t. fi{x) < for \x\ < 6, but fi is not always negative; 
Dl) there exists c> s.t.: 

Fi{x) > —c if X > and Fi{x) < c ifx<0; 

D2) limsup3,^±^ [G{x) ± Fi{x)] = +oo; 
D3) lima;^±oo Jq /2(s) ds = l±> -oo, 
E) there exists A' > s.t. either: 
El) f2ix) >Oforx< -A', 



E2) 
or 



lim sup2,^_^ 



Fiix) + 



El') f2{x) <Oforx>A', 



E2') limsup^ 



F,ix)- 



L_ < +CX); 



= Li > —oo. 



Then system ((1.3|l has at least a periodic orbit. 

Certainly this result can be somehow, directly obtained after transforming the 
initial equation into a Lienard's one, but however we point out that in our for- 
mulation based on a geometric study of the phase-plane, assumptions keep their 
geometrical meaning, as has been used in |llj. Moreover in some applications, as 
for instance Theorem 13. 31 they can be easily verified. 

Some examples of limit cycles for the general case, when there is not global 
boundedness of solutions, are also presented. 

In Section 0] we consider the problem of uniqueness of limit cycles. A classical 
result due to Massera is generalized as to include the following equation: 



(1.5) 



N 

E 

/=0 



+ V/2;+i(^)i"+'+a^ = 0, 



for some positive integer N . The Classical Massera result is recovered when iV = 0. 

Theorem 1.3 (Existence and Uniqueness). Let f2i+i{x), I G {0, ...,iV}, N > 

1, verify assumption A and moreover: 

LI) there exists 5 > s.t. fi{x) < for \x\ < 5; 
L2) f2i+i{x) > 0, ; e {1, . . . ,iV}, for all x; 

L3) f2i+i{x), I G {0,...,Af}, is monotone increasing (respectively decreasing) 
for a: > (respectively x <Q). 
Then (|1.5|l has a unique globally attracting limit cycle. 

As far as we know there are no other generalizations of classical Massera's result 
and the problem of adapting his geometrical ideas to a more general situation, still 
remain an open question. For this reason we believe that this particular result is 
in some way significant. 



4 



TIMOTEO CARLETTI, LILIA ROSATI AND GABRIELE VILLARI 



Remark 1.4. Assumptions LI and L3 for fi{x) don't ensure the existence of the 
limit cycle in the Massera's case (N ~ Q); on the contrary our assumptions guar- 
antee that (|1.5|l always has a limit cycle. For instance if fi{x) = (x^ — l)e~^ 
in Massera's case the limit cycle doesn't exist ^ but adding a "perturbation", small 
as we want, say f3{x) = tx^ jix? + 1) with e > 0, we ensure existence, and thus 
uniqueness of the limit cycle. 

2. A NON-EXISTENCE RESULT 

The aim of this section is to prove the non-existence result Theorem ll.ll Through- 
out this section regularity hypotheses to guarantee existence and uniqueness of the 
Cauchy initial problem, will be assumed; also sign assumption B will hold. 

We will prove this result by considering firstly the simplest case of Hl.r>|l dealing 
with only two terms, and secondly the general situation involving several terms. 

In the phase-plane equation H1.3|) can be rewritten as a first order differential 
system: 



u = V 

,2 



^ = -fi{u)v - f2{u)v^ - g{u) , 

thus using the well known transformation 



(2.2) 



X — u 



\y = - /o" h{s)exp\^ /2(r) drj ds^ exp\^ - /g" /2(r) drj 
with the rescaling of time: 

(2.3) dr = exp\^ - f /2(s) dsj , 



dt, 

system H2.1II can be brought into a Lienard one (still denoting by i = dx/dr) 
(2.4) 



X = y — F{x) 
y = -9{x) , 



where: 

F{x) = J fi{s)exp^ J /2(r)dr|ds and g{x) ^ g{x)exp^2 J f2{r)dr^. 

A large number of results are based on determine hypotheses on F and g to ensure 
existence, uniqueness or non-existence of limit cycles for system (|2.4() . and then 
to transport them to the original system. We will adopt a different point of view, 
considering system H2.1|l as a "perturbation" of the system: 

{ii = V 
V ^ - f2{u)v'^ - g{u) , 

and exploiting the symmetry of the last one w.r.t. v ^ —v. 



^This claim follows easily studying the Massera case in the associate Lienard plane, where the 
flow through circles {x^ + = r^} is —x fi (s) ds, which doesn't change sign with this choice 
of 
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The "unperturbed" system, i.e. with fi{x) 
tonian one: 



0, is transformed into an Hamil- 



(2.5) 



y 



whose Hamilton function is: H{x,y) = + G{x), where G{x) = Jq g{s) ds. We 
remark that the origin is a local center. It can be a global one according to the 
behavior of G{x) for large |a;|. 

Remark 2.1. Restricting to the class of polynomials the origin is a global center 
if and only if /2 (x) is polynomial of odd degree with positive leading coefficient. In 
all the remaining cases the origin is a local center and at least a separatrix appears. 

The transformation 1)2.2(1 is invertible, hence one can bring back to the original 
plane the energy curves {H = const}: they will keep the same topology being only 
slightly distorted in the phase-plane (see for example Figure QJ. 




Figure 1. fi{x) = + 1, f2{x) = -x^, g{x) 
{H = A} with A e {0.1, 0.3, 0.5, 0.7, 0.9, 1.1}. 



X, energy levels 



Let us now suppose the perturbation fi{x) ^ 0. We will show, using simple com- 
parison arguments, that system l|2.1|) cannot have limit cycles if /i never changes 
sign. Let us suppose fi{x) > for all x, and let us compare the slopes of sys- 
tem 1(2. l|l . when /i > 0, but not identically zero, say case A, and fi=0 say case 



B: 



dv 
du 



dv 



< 



dv 
du 



A du 

Hence trajectories in case A enter trajectories of case B and no limit cycle can exist. 
This concludes the proof for 

We point out that this situation occurs both if 1(2. 5|l has a global or local center, 
namely this situation is not the perturbation of a global center. 
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We remark that this is not a new result, in fact if /i keep constant sign, then 
F vanishes only at x = 0, and it is positive (respectively negative) for x > 
(respectively x < 0), thus Lienard system l|2.4|l has not limit cycle and so does the 
initial system. But our generalization, obtained from Remark 12 . 21 and Remark 12. 31 
is to the best of our knowledge new, and it cannot be obtained "passing through 
some Lienard system" because the previous transformation cannot be anymore 
extended. 



Remark 2.2. We can allow terms f2i{x)x^^ into (|1.3|) 
(2.6) x + fi{x)x 



M 

J2f2i{x)x'' +g{x) 
1=1 



0. 



For such kind of systems there is not generically a transformation mapping them 
into some Lienard systems. Nevertheless we can prove our non-existence result. 
A sketch of the proof of this claim is as follows: assume fi = and consider 
systems H2.6|l in the phase-plane. Then there is a trivial symmetry y i— > —y which 
allows to prove the existence of a local center. Thus, once we add e perturbation 
fi{x) > we can use the same reasoning as before, comparing slopes, to prove the 
non-existence of limit cycles. 

The following remark concludes the proof of Theorem ll.il 

Remark 2.3. Provided f2i-i{x) > 0, Z G {1, . . . ,iV} for all x, we can generalize 
further (|2.6|l by considering: 



N 



M 

i + Y. hi-i{x)±^''^ + f2i{x)x^' + g{x) = 



1=1 



1=1 



The proof is straightforward and we omit it. 

3. Existence of limit cycles 

In this part we will prove two existence results; the first one, Theorem 11.21 is 
more general for the involved functions but deals with system H1.3|l . the second 
one. Theorem 13.31 holds only for polynomials but for a slightly more general equa- 
tion (ESI). 

Theorem ll.2l fExistencel. Assume hypotheses A and B to hold, let Fi{x) = 
lo fi{s)ds , G{x) = /(^ g{s)ds and assume: 

C) there exists 5 > Q s.t. fi{x) < for \x\ < 6, but fi is not always negative; 
Dl) there exists c> s.t.: 



Fi{x) > —c if X > and 
D2) limsup^^i^ [G{x) ± Fi{x)] = +oo; 
D3) lim^^ioo Jq /2(s) ds = l±> -oo, 
E) there exists A' > s.t. either: 
El) f2{x) >Oforx< -A', 



Fi{x) < c ifx<0; 



E2) limsup,^_^ [F.ix) + 
or 

El') f2ix) <Oforx>A', 



E2') limsup,^+^ [f,{x) - 



L_ < +oo; 



= L+ > —oo. 
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Then system ()1.3(l has at least a periodic orbit. 

Remark 3.1 (About hypothesis D3). We point out that hypothesis D3 can he 
removed, but in that case some growth conditions on fi{x) and/or g{x) have to 
be assumed. The exact ones would be clear in the proof of Theorem M.SX Roughly 
speaking if f2{s)ds diverges to — oo as some odd power, x^^^^ , for x — > —oo, 

then fi{x) has to grow faster than e^ ^ to "compensate" this divergence. If not 
the intersection property with the oo-isocline is not guarantee, thus there can be 
non-winding trajectories. 

Remark 3.2 (Polynomial case). The sign assumptions on fi (conditions El or 
El') and hypothesis D3 cannot he simultaneously verified if f2(x) is a polynomial. 
Hence Theorem \l.'2\ cannot be applied to the important class of polynomials. Nev- 
ertheless, we will provide in § a result ( Theorem for the following second 
order differential equation with polynomial coefficients: 

p{x)x + p{x)qx{x)x + q2{x)x^ + r[x) = , 

which is strictly related to (|1.3|l . 

3.1. Proof of Theorem II. 2L As already remarked in the previous section, sys- 
tem (|1.3|) rewrites, in the phase-plane: 



(3.1 



V 

,2 



\v = -fi{u)v ~ f2{u)v^ - g{u) . 

Let us introduce a new systems of coordinates x — u and y — v + Fi {u) , such 
that the previous phase-plane system rewrites as: 

J i = y - Fi{x) 



(3.2) 



-f2{x){y-F,{x)f~g{x). 



Hypothesis C ensures that the origin is a repellor, in fact looking at the flow 
through the ovals Or = + G{x) = r^} and r is a small positive number, we 

get: 

^Or ^ = ~9{x)F,{x) - yf2{x) {y - F,{x) f , 
at flow 

and the claim follows remarking that the last term in the right hand side is of higher 
order than —g{x)Fi{x) close enough to the origin, thus the sign is determined by 
the first one. 

Thus to prove the existence of a closed trajectory is enough to prove the existence 
of a solution spiraling toward the origin. 

First of all we want to guarantee that all trajectories intersect the curve y = 
Fi (x) . For Lienard system one has a necessary and sufficient condition to guarantee 
such property, moreover as already remarked in the previous section, system (I3.2|l 
is equivalent, through the transformation 1)2. 2|l and the rescaling of time (|2.3I) . to a 
Lienard one : 



(3.3) 
where: 

m 



fi{s)exp^j f2{r)dr^ds and g{€) = 9{^)(^xp^2 j /2(r)dr| 
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We observe that the curve y — Fi{x) corresponds in the phase-plane (|3.HI to the 
curve V = 0, which in the Lienard plane 1)3. 3|l transform into rj — -F(C), hence 
intersection with y = Fi{x) occurs if and only if intersection with rj — F{£^) occurs. 

Let G(0 =J^g{s)ds. Ass uming conditions B and Dl one can apply the follow- 
ing result jTU) : 

• for any > and 770 > F{£^q), the solution of H3.3|l passing through (Co, ??o) 
will intersect rj = F{^) at some ^f', F{^')j , with ^' > ^o, if and only if: 

(3.4) limsup [g(C) +/'(0 

• for any Co < and 770 < F{£.o), the solution of 1)3. 3|l passing through (Co, ?7o) 
will intersect rj = F{^) at some ^C', ^(O) , with C' < Co, if and only if: 

(3.5) limsup [g(C) -f'(C) 

We claim that from hypotheses D2 and D3 conditions (|3.4|l and (|3.5|l follow. Let 
us look for example to G(C) assuming /q^°° /2(s) ds — 1+ to be finite and positive, 
then by definition of limit there exists M > s.t.: 

for aU x>M. Let C > M, by definition: 

G(C)=/ g{x)e-^^S h{s)ds^^^ 5(x)e-2/o"/2(-)'i- dx, 

thus the first term is a constant, depending on M, whereas the second can be 
estimated by: 



g{x)dx< / 5(x)e-2/o"-^2(s)ds^2; < e-'+ / g{x)dx. 
M Jm Jm 



Finally putting everything together we obtain, for all C > M, the bound for G(C): 
G(M) + e-3'+ (G(C) - G(M)) < G(C) < G(A/) + e-'+ (G(C) - G(M )) , 

which implies that G(C) and G(C) have the same behavior for large C- 
The other cases can be handled similarly and we omit them. 
In this way we have the intersection property: trajectory of (|3.2|l passing through 

(a;o,yo), where a:o > and yo > Fi{xq), will intersect the curve y = Fi{x) at some 

{x' , Fi{x')) with x' > xq. And similarly for the one starting at (a;o,yo), where 

xo < and yo < Fi{xo). 

To conclude the proof we must control that not all trajectories escape. To do 

this we will study the zero-isocline curves which are explicitly given by: 



9{x) 



Let us assume hypothesis El to hold (the case for El' can be handle similarly 
and we will omit it). Because of the sign assumptions on 72(2;) and g{x), y^ocix) 
are well defined for x belonging to (~oo, —A') U X, where T is the union of open 
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disjoint intervals, whose endpoints are zeros ^ of f2{x). This gives rise to (possible 
several) branches for the zero-isocline curves. Let y^j:^{x) be the leftmost one above 
y = Fi{x), defined for x < —A'. Such branch can be used, assuming hypothesis 
E2, to prove the existence of a spiraling orbit toward the origin. 

Before to prove this fact, we have to control the continuability of solutions. We 
claim that for all a < /3 < (respectively < a' < (3') and for all y > (respectively 
y' < 0), the solution of (|3.2|l with initial datum (a, y) will intersect the line x = (3 in 
finite time (respectively the solution of (|3.2(l with initial datum (/?', y') will intersect 
the line x = a'). To see this just consider the system in the phase-plane and recall 
that in this transformation vertical lines are mapped on vertical lines. The slope of 
such a solution is dv/du = —fi{u) — f2{u)v — g{u)/v, hence if is large enough, 
this slope is sub-linear, giving rise to the continuability from a to (3. Moreover, the 
explicit solution is bounded by some exponential function which is the responsible 
for the "strong" distortion of the limit cycles shape (see Figure OJ. 

We are now able to conclude our proof. Let us consider an orbit starting at 
xq < —A', 2/0 = Viefti^o), recall that we are assuming conditions El and E2 
to hold, if instead of those we were assuming hypotheses El' and E2' we should 
start from the rightmost branch of the zero-isocline below y = Fi (x) from some 
a;o > A'. Such trajectory will intersect the y axis by the previous discussion about 
continuability (that's arc AB in Figure EJ. Then, thanks to condition H3.4(l it will 
reach the curve y = Fi{x) passing eventually through branches of zero-isocline 
(arcs BC and CD in Figure[2Jl- On y = Fi{x), x changes sign and the solution 
goes back till intersecting again the y axis at some negative ordinate (arcs DE 
and EF in Figure 01 • Using now (|3.5|l we are able to prove that this trajectory 
will reach again the curve y = Fi{x) at some a; < 0, where again x changes sign. 
Before reaching this curve, it has to intersect a branch of zero-isocline below the 
curve y = Fi{x), hence y > 0. Because this trajectory cannot intersect the leftmost 
branch of y'^^j^{x) upon which the vector field points to the right, it will reach anew 
the quadrant from which it started spiraling toward the origin when doing a tour 
around it. Thus it can be used as outer boundary of a Poincare-Bendixon domain 
and this concludes the proof. 

3.2. The polynomial case. As already remarked Theorem 11.21 cannot be applied 
to the class of polynomials, in fact the sign assumptions (conditions El or El') and 
hypothesis D3 cannot be simultaneously verified by polynomials. Nevertheless we 
can prove a similar result for a slightly different second order differential equation 
with polynomial coefficients. 

Theorem 3.3. The second order differential equation: 

(3.6) p{x)x + p{x)qi{x)x + q2{x)x^ + r{x) = , 

where p, {q)j=i^2 o-nd r are real polynomials of the x variable, has at least a limit 
cycle provided the following conditions hold: 

HI) p{x) > for all x, qi{0) < and xr{x) > for x ^ Q; 

H2) qi{x) has even degree and the leading coefficient is positive; 

H3) degpix) >degq2{x)+2; 

^If /2(0) = 0, because also g(0) = 0, the point will be an endpoint of some interval in X 
according to the behavior of g{x) / f2{x) close to x = 0. 
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Figure 2. The zero-isocline curves and the oo-isochne curve y 
Fi[x), with a trajectory spirahng toward the origin. 



H4) if q2{x) has odd degree then the leading coefficient must be negative; 
H5) degr(a;) < 2degqi{x) + degq2{x) + 1. 

We observe that the special form for the coefficient of x is not a restrictive one, 
we choose it only to state hypotheses of Theorem 13.31 in a simplest and more clear 
way than using some polynomial qi{x) instead of the factorization p{x)qi{x). In 
this way, the applicability of the present result is straightforward: just look at the 
polynomials and compare their degrees. 

Proof. Because p{x) never vanishes we can divide by it, then we set fi(x) = 
qi{x), f2{x) = and g{x) = ^[fy- The proof will be complete once we show 

that /i, /2 and g verify hypotheses of Theorem ll.2l 

Existence and uniqueness of the Cauchy initial value problem associates to (|3.6|) 
is trivially guaranteed because we are dealing with polynomials. 

Signs assumption HI ensures that hypothesis B holds and by continuity of qi (x) 
also the first part of assumption C is verified. Also the second part of it holds, 
because qi has even degree and positive leading coefficient. 

Hypotheses Dl and D2 are verified because the primitive of qi{x) is an odd 
degree polynomial with positive leading coefficient and moreover the primitive of 
g{x) vanishing at x = is not negative. 

By H3 it follows that the integral q2{s)/p{s) ds is well defined and finite for 
all X. 

Let a (respectively /3) be the smallest (respectively largest) zero of q2{x). If q2 
has even degree and positive leading coefficient, then q2{x), hence f2{x), is positive 
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for all X < a. Moreover let Qi{x) — qi{s) ds, then for x — > — oo: 



(^) _ ^dcg 91 +1 L ^ . /^I ^1 ^"2'°' -dog 91-1 



where a, b and c are the positive ^ leading coefficient of respectively Qi, r and q2. 
Hence assuming H5 hypothesis E2 is verified. 

If (72 has even degree and negative leading coefficient, then q2{x), hence f2ix), is 
negative for all x > (3. Moreover for x — > +00 one has: 



and again H5 ensures hypothesis E2' to hold. 

Finally if (72 has odd degree and negative leading coefficient, then both couples 
of hypotheses El, E2 and El', E2' hold with A' = max{|Q!|, |/3|}. □ 

In the following Figure Owe show some "exotic" limit cycle whose existence is 
guaranteed by this last Theorem. 



1 y\ \ 
1 ^\ \ 

/^\' 

/^\\ 
I 

• No 



\ \--/ / 




/ / I 



Figure 3. We show an exotic limit cycle, whose existence is 
proved by Theorem 13.31 -pix) = (a;/3)^ + 1, q\{x) = x"^ — 1, 
92(2;) — X* — Ax'^ + 1 and r{x) — x. 



To conclude this section let us discuss briefly the general situation where there are 
few results about existence of limit cycles. This is because, as already mentioned, 
one cannot always have continuability of solutions, hence the Poincare-Bendixon 
Theorem, as usually done in Literature, cannot be applied "in large" . Nevertheless, 

■^Assumption on the sign of r implies tliat r cannot have even degree nor negative leading 
coefficient. 
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using standard Hopf's bifurcation argument, we are able to produce examples with 
small amplitude limit cycles, even if global continuability is not allowed. 
To fix the ideas, let us consider the following example: 

X + {ax^ - b)x + (a;^ + l)x^ + x^x^ + .t = ; 

phase-space investigation shows that the slope is superlinear, hence not all trajec- 
tories are continuable. Near the origin the dynamics is dominated by the x term, if 
a > and 6 = the origin is a local attractor. For positive, suitable small h a Hopf's 
bifurcation takes place: the origin changes its stability and a small amplitude limit 
cycle appears, while "in the large" orbits behavior doesn't change too much. 

The idea can be generalized as to include other situations, where always small 
amplitude limit cycles exist and there is not global continuability of solutions. 



4. A GENERALIZATION OF A THEOREM BY MASSERA 

4.1. Introduction. Massera's Theorem (see but also |^ pages 307-308) con- 
cerns the problem of existence and uniqueness of limit cycles for the second order 
differential equation (|1.5() . which in phase-space can be rewritten as: 



(4.1) 



X =y 

y =-fi{x)y-x. 



Under the assumptions: fi{x) < for \x\ < 6, for some positive (5, and /i increasing 
(respectively decreasing) for a; > (respectively a; < 0), Massera proved that if a 
limit cycle exists, then it is globally attracting, hence unique. The proof is based on 
the "geometrical" construction of the so called geodesic system, according to Z.-F. 
Zhang 1121. 

The previous assumptions on /i (x) do not guarantee the existence of the limit 
cycle; an easy calculation (see for instance footnote H at page^J shows that sys- 
tem H4.1|l has not limit cycles at all, if for all 6i < S2 one has: J^^ fi{x) dx < 0. 
We will see in a while that this situation cannot occur in our case. See Remark ll. 41 
and Figure 13 for an explicit example. 

Generalizations of this Theorem, for instance by replacing y = —fi{x)y — x 
with y — —fi{x)y — g{x), are not available. Here we propose a new point of view 
considering the following system: 



(4.2) 



= y 

^~h{x)y~j:lif2i+iix)y''+' 



under the additional hypothesis: f2i+i{x) > for all x and f2i+i{x) is increasing 
(respectively decreasing) for x > (respectively x < 0). 

We point out that system 14.2|l admits a limit cycle even if system H4.1|l doesn't, 
moreover if the classical Massera system has a unique limit cycle then it must 
contain in its interior the one of 1)4. 2|l . 



4.2. Existence of limit cycles. Let us consider the circles Cr = {{x,y) S IR^ : 
x'^ + y'^ — r^}, we claim that the flow of (|4.2(l is transversal (more precisely it points 
outward) to Cr for r small enough, hence the origin is a repellor. 
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To prove this claim let us consider the angle between Cr and the vector field 
X = {y, -fi{x)y - /2i+i(a;)y^'+^ - x) valued at C^: 



= < VCr,X > 



N 

1=1 



thus if r > is small enough, fi is negative and the remaining positive terms cannot 
compensate it. 

On the other hand if r is large enough then ar < 0, in fact its leading term, 
— 2y^^+^/2Ar+i(a:), is negative because of the positiveness of f2N+i{x). Hence by 
Poincare-Bendixon's Theorem we have constructed the inner and outer boundaries 
of an invariant domain which must contain at least one limit cycle. 



4.3. The cycle is star— shaped. A fundamental assumption for the proof of 
Massera's Theorem is that the limit cycle must be star-shaped w.r.t. the origin. 
Z.-F. Zhang (see T? pages 246-253) pointed out, that generalizations of Massera's 
Theorem following the geometrical ideas of geodetic systems, must rely on this 
assumptions. 

The aim of this section is to prove that this condition is verified for system (|4.2|l . 

Let assume the limit cycle to be not star-shaped w.r.t. the origin, thus there 
is at least a half-line from the origin intersecting the cycle at least at two points. 
Starting from the point A on the y-axis the vector field points right-down, i > 
and 2/ < if a; is small enough, hence the geometry must be as in Figure^ 




Figure 4. Part of a non star-shaped limit cycle. 



Let A — PO/BO and let us construct the arc RP homothetic to AB^ with 
ratio A. Comparing the slopes of the vector field on these two arcs we will show a 
contradiction, hence the limit cycle must be star-shaped. Take any point {xi^yi) 
on the arc AB and its homothetic one {x2,y2)j let y = mx be the straight line 
passing through these two points. We can assume X2 > xi > 0. Let us compare 
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the slopes at {xi,mxi) and {x2,mx2)'- 

— {X2,mx2) = -fi{x2) -^f2i+i{x2){mx2f'^ - m 

1=1 

(4.3) < -fijxi) - f2i+i{xi)(mxi)'^^ - m = —{xi,mxi) , 

1=1 

where we used the fact that the functions x ^ — (ma;)^'/2i+i(x), I £ {0, . . . , iV} 
are decreasing for x > Q. This imphes that orbits starting on the arc RP move 
toward the arc AB. But this gives rise to a contradiction, in fact for these orbits 
we have a; > 0, they cannot accumulate to a point somewhere in the generalized 
quadrilateral with vertices S, P and (2:2, j/2) because the origin is the only 

fixed point, thus they must intersect the limit cycle somewhere between points 
A and P, which is impossible because the system is autonomous and Cauchy's 
uniqueness and existence assumptions are guaranteed. 




Figure 5. The hmit cycle for system (|4.2ll with: fi[x) — (x^ - 
l)e^^ , /3(a;) — x^/50(a;^ + 1) (solid line), and the spiraling orbit 
for the classical Massera H4.1|l with initial datum (0,0.01) (dashed 
line). 



4.4. Uniqueness of limit cycle. In this paragraph we will prove uniqueness of 
the limit cycle, by showing its hyperbolicity. In our proof we will follow closely the 
original ideas of Massera. 

Let us take a limit cycle, 71, and consider its homothetic by a factor /i > 1. 
Take a point (xi, y\) G 71 and its homothetic one (x2, 2/2) G 72, let y = rax be the 
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straight line passing through these two points. Comparing again the slopes, as we 
did in H4.3(l . we get: 

dy , \ ^ dy 

dx ax 

for all xi and X2- Here the equal sign appears because we are considering also 
xi = 0. Hence trajectories starting on 72 will go inward and accumulating to 71. 

Similarly, considering 73 homothetic of 71 by a factor < < 1, and comparing 
again slopes we will obtain that 71 is internally attracting. This proves that 71 is 
globally hyperbolic, hence unique. 
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